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SYSTEM AND METHOD FOR CUMULANT-BASED GEOLOCATION OF 
COOPERATIVE AND NON-COOPERATIVE RF TRANSMITTERS 

Related Applications 

[0001] The present application is related to and co-pending with commonly-assigned 
U.S. Patent Application S.N. 10/360,631 entitled "Blind Source Separation Utilizing A Spatial 
Fourth Order Cumulant Matrix Pencil", filed on 10 February 2003, the disclosure of which is 
hereby incorporated herein by reference. 

[0002] The present application is related to and co-pending with commonly-assigned 
U.S. Patent Application S.N. 10/400,486 entitled "Method And System For Waveform 
Independent Covert Commimications", filed on 28 March 2003, the disclosure of which is 
hereby incorporated herein by reference. 

[0003] The present application is co-pending with and claims benefit of U.S. Provisional 
Patent Application Serial No. 60/458,038 entitled "Cooperative SIGINT for Covert 
Communication and Location Provisional", filed on 28 March 2003, the entirety of which is 
hereby incorporated herein by reference. 

Government License Rights 
[0004] The U.S. government has a paid-up license in this invention and the right in 
limited circumstances to require the patent owner to license others on reasonable terms as 
provided for by the terms of Contract No. NRO000-02-C-0389 awarded by the National 
Reconnaissance Office. 

Background 

[0005] Information regarding location of a source such as for surveillance or combat 
search and rescue can be degraded in value if detected by unfriendly entities, such as enemy 
forces in the case of a downed pilot or a marked terrorist under surveillance 



WSH\96179.1 



1 



GCSD-1448 



HAR62-013 



[0006] Intentional detection of the signal or message can be accomplished in military 
systems that use specially designed electronic support measures (ESM) receivers. These ESM 
receivers are often found in signal intelligence (SIGINT) applications. In commercial 
applications, devices employed by service providers (i.e. spectral monitors, error rate testers) can 
be used to detect intrusion on their spectral allocation. Inadvertent detection can also occur, such 
as when a user or service provider notices degradation in link performance (e.g., video quality, 
audio quality, or increased bit error rate). 

[0007] The term covert also implies the additional goals of evading interception and 
exploitation by unintended receivers. Interception is the measurement of waveform features or 
parameters useful for classifying/identifying a transmitter and/or the waveform type and/or 
deriving information useful for denying (i.e. jamming) the communication. Exploitation is 
processing a signal by an unintended receiver in the attempt to locate the transmitter and/or 
recover the message content. In the broad literature on covert communications these 
characteristics as applied to transmitted information signals are referred to as low probability of 
detection (LPD), low probability of intercept (LPI), and/or low probability of exploitation (LPE) 
by an unintended receiver. 

[0008] Given the desirability to transmit messages covertly, it is helpful to understand 
considerations that enhance or degrade LPD, LPI and LPE. An unintended receiver such as the 
receiver 103 in Figure 1 with the goal of detecting a covert communication must reliably 
differentiate between the binary noise-only and signal-plus-noise hypothesis. As is known to 
those of skill in the art, for an unintended receiver the signal detection process is typically based 
on an energy threshold. The energy the receiver measures is given by E^^^ = Pavg^xmit • Where 
under general conditions the power Pavg is the received covert signal power S plus intemal 
receiver noise power N. Hence, E^^^ = (^S + N)T^^i^ . If the signal power used to communicate 

is only a small fraction of the receiver noise, S«N, it is extremely difficult for the unintended 
receiver to reliably detect the presence of the covert signal because the total energy detected will 
only be marginally greater than the noise-only (5=0) case. 

[0009] Minimizing transmit power has two direct system benefits. First, the total signal 
power used will be a small fraction of the total noise power present in the same band. Thus, if 
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the message is limited in time duration, the total energy measured by an unintended receiver 106, 
which may be an ESM receiver, is indistinguishable from a noise-only environment. Since ESM 
receivers are often of energy threshold type, there is an obvious trade-off of average power for 
time duration in order for a signal to remain imdetectable. Second, the low transmit power 
scenario enables usage by transmitters with very limited power supplies (i.e. batteries). 

[00101 Therefore, as naturally arise in military environments such as depicted in 
Figure 2, there is a need for a low power message system and method, covert or otherwise, such 
as covert communications for Intel or Special Forces, "stealth" IFF for low observable ground 
vehicles, and combat search and rescue (CSAR). There is also such a need in a nimiber of 
civilian or public safety applications as well, such as asset tracking/location or "lost child" 
detection/location and surveillance. In particular in these latter-described applications it may be 
particularly desirable to receive both a message and location the source of the message. 

[0011] As mentioned above it is often of interest to geolocate signals, particularly those 
that may de designed for LPI/LPS. These included spread spectrum signals, spread spectrum 
signals are intentionally low power as previously discussed, and these signal can also be co- 
channel with many other signals of similar type, which makes geolocation by prior art methods 
and systems ineffective. 

[0012] Embodiments of the present inventive system and method address the above 
needs while requiring only an extremely low power signal. The geolocation needs are 
specifically addressed by estimating two cxunulant matrices, and performing generalized 
eigenvalue decomposition (GEVD) of the resulting matrix pencil. The GEVD provides 
eigenvectors orthogonal to the incoming steering vectors, save one. Exploiting this property 
allows estimating of steering vectors for each incoming signal. From the steering vectors it is 
easy to arrive at AOA or geolocation. The embodiments enable geolocation signals that are 
below thermal noise and in co-channel environments. 

[0013] These and other advantages of the disclosed subject matter will be readily 
apparent to one skilled in the art to which the disclosure pertains from a perusal or the claims, the 
appended drawings, and the following detailed description of the preferred embodiments. 
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Brief Description of the Drawings 
[0014] Figure 1 is a general representation of a prior art approach to LPD. 
[0015] Figure 2 depiction of covert scenarios. 

[0016] Figure 3 is a representation of an embodiment of a waveform independent covert 
communication system. 

[0017] Figure 4 is a depiction of a binary symbol message 101011 according to an 
embodiment of the invention. 

[0018] Figure 5 a is a representation of a spectrogram of a transmission of the binary 
symbol message depicted in Figure 4 and revealing the symbol boundaries of the message 
according to an embodiment of the invention. 

[0019] Figure 5b is a representation of a spectrogram of a transmission of the binary 
symbol message depicted in Figure 4 as might been seem by a typical ESM/SIGINT spectral 
monitor. 

[0020] Figure 6a is a generalized schematic representation of embodiments of the 
invention. 

[0021] Figure 6b is a schematic representation of an embodiment of the invention with 
dual signal sources. 

[0022] Figure 6c is a schematic representation of an embodiment of the invention with 
dual temporal filters. 

[0023] Figure 6d is a schematic representation of an embodiment of the invention with 
dual non-Gaussian noise generators. 

[0024] Figure 6e is a schematic representation of an embodiment of the invention with 
dual non-Gaussian noise generator and temporal filter branches. 
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[0025] Figure 6f is a schematic representation of an embodiment of the invention with 
dual temporal filters and dual signal sources. 

[0026] Figure 7 is a schematic representation of a Laplacian noise generator with 
multiple signal sources according to an embodiment of the invention. 

[0027] Figure 8 is a representation of eigenvalue tracks for a firequency shifter laplacian 
noise waveform for message "101011" according to an embodiment of the invention. 

[0028] Figure 9 is a representation of bloclc-to-block eigenvalue correlations. 

[0029] Figure 10 is a schematic representation of a message recovery system with spatial 
information according to an embodiment of the invention. 

[0030] Figure 1 1 is a flow diagram for covert communication via message recovery with 
spatial information according to an embodiment of the invention. 

[0031] Figure 12 is a schematic representation of a message recovery system without 
spatial information according to an embodiment of the invention. 

[0032] Figure 13 is a flow diagram for covert communication via message recovery 
without spatial information according to an embodiment of the invention. 

[0033] Figure 14 is a flow diagram for covert commimication via encoded eigenvalues 
according to an embodiment of the invention. 

[0034] Figure 15 is a representation of a binary message sequence 101011 encoded using 
different carrier waveforms. 

[0035] Figure 16 is a representative diagram for geo-locating a transmitter using higher 
order cimiulants. 

[0036] Figure 17 is a flow diagram for tracking a cooperative transmitter. 
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Detailed Description 

[0037] A useful feature of embodiments described herein is the ability to geolocate a 
remote transmitter independent of the waveform or protocol used. The individual waveforms 
used to communicate are in a sense superfluous or independent of the message to be conveyed. 
This is a significant advantage for ubiquitous application, allowing for parasitic use of present 
communication infrastructure and devices. Thus there are few restrictions on the pairing 
between potential covert transmitters and the intended receiver using the disclosed covert 
communication methods and apparatus because of the independence of the information transfer 
on the "carrier waveform". This is unlike prior art systems where the receivers designed or 
instantiated for a certain signal type cannot accurately recover the message if the receiver is 
presented with another signal type. However, embodiments of the present invention by contrast 
can function equally well for any waveform, and the location system does not require any a 
priori knowledge of the "carrier waveform". In fact, embodiments of the covert transmitter can 
be waveform agile without informing the intended receiver. 

[0038] The embodiments herein are predicated on selecting and transmitting carrier 
waveforms with unique higher order spatial statistics. Such higher order statistics include 2"^ 
order spatial moments and 4* order spatial cumulants. . The primary restriction is that the 
receiver and transmitter must use the same "codebook" of time durations and alphabet. . 

[0039] To recover the message information, the waveform is not conventionally 
demodulated. Rather, a straightforward block or batch estimation algorithm estimates the 
generalized eigenvalues of the SFOCMP for each signal in the receiver field-of-view (FOV) over 
time using only the array output. For this discussion we assume the data has been digitized 
appropriately. The sizing for the block processing (e.g., the block of contiguous array 
observations, sometimes known as "snapshots") is dependent on several factors. Chiefly we 
must ensure that each block has enough sample support so that the eigenvalue estimates fi-om the 
GEVD of the SFOCMP in each block over a symbol duration nominally match. This means that 
estimation error is negligible. Accordingly, changes in the eigenstructure can be reliably 
detected, and this change indicates a symbol boundary. The degree to which a nominal match is 
required within a block depends on the complexity of the signal environment (e.g., extraneous 
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co-channel signals), the communication errors (e.g., partially received messages) tolerable in a 
given application, and the receiver processing resources to recover the message in a timely 
manner. 

[0040] Further, in practical situations, as power sources become impaired (e.g., batteries 
running low on power), transmitted waveforms become increasingly distorted. This situation 
limits the effectiveness of matched filtering as used in prior art systems, since the concept of 
matched filter relies on knowledge of the transmitted waveform in the receiver. Embodiments of 
the inventive technique are impervious to such distortion since it is the duration and not the 
actual value of the eigenvalues of the SFOCMP that matter. So as long as the eigenstructure 
characteristic of the distorted signals is nominally constant during a message symbol, the 
inventive system and method is robust as to degraded transmitter performance. Therefore, the 
present inventive system and method will operate successfiiUy imder conditions that would 
normally be detrimental to conventional systems. The use of lower order matrix pencils are also 
contemplated by the present inventive system and method. 

[0041] The digital data can also include fi*aming or formatting of the message. Typical 
of the fi-aming would be start/stop and data fields. Though other fields can be used as needs 
dictate. This fi-aming structure can ensure that the receiver can reliably find the beginning of a 
message for synchronization. Source coding or compression could also be applied to the 
incoming data stream to reduce the required bandwidth. The user may also encrypt the data to 
protect it by an optional encryption device of suitable complexity. This data is then output to the 
forward error correction (FEC) module, which currently is envisioned as applying block coding. 
The coding is usefiil to aid the receiver in resolving message ambiguities say caused by fades or 
other unresolved time coincident measurements, which in this system could be processed as 
"erasures" up to the correction capability of the code Thus, the potential for a message protocol 
using automated repetition of the message might be advantageous as error patterns in each 
transmission will likely be different. 

[0042] The "carrier waveform" has its fourth-order characteristic modified according to 
the control by the M-ary alphabet. There is no limit to the strategies potentially adopted by the 
covert transmitter for this operation, so long as the characteristic is measurable reliably by the 
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receiver and it conforms to the time duration and alphabet size assumed by the receiver. To this 
end, there will be some minimum duration and maximum duration for a symbol, and a preferred 
duration increment for each symbol in the alphabet. The exact choices depend on application, 
however, making the durations too disparate can negatively affect data rate limiting this 
technique to lower data rate applications. It is desirable to provide durations that are easy to 
resolve into the M-ary symbols. 

[0043] The minimum duration and duration increment must be such that synchronizing 
the data block boundaries used in the receiver to that of the symbol timing in the transmitter 302 
is not relied upon. It is desirable that the covert transmitter 302 use a fundamental signaling 
period of several (e.g., 10) "receiver block" durations for the minimum signal, and may have a 
signaling duration increment of the same size to define the alphabet. However other choices are 
applicable depending on the particular implementation and application. The goal is to make the 
time duration alphabet as disparate as possible while meeting performance objectives (e.g., data 
throughput). Sample data rate computations can be determined as shown below. 

[0044] Defining S as the array snapshots/block, "b" blocks for the minimum length 
symbol, "B" blocks for the maximum length symbol, and "R" as ADC (analog to digital 
converter) conversion rate in the receiver, the minimum and maximum symbol durations for a 
binary alphabet are: 

[0045] T,y„(min) = S*b*l/R = 5,000*10*10-^ 50 x lO^sec 

[0046] TUmax)= S*B*1/R = 5,000*20*10-9= 100 x lO^sec 

[0047] The values b=10 and 3=20 along with S=5,000 and R= 1 Gsamples/sec are 
subject to implementation choice, and used here for illustration only. Assuming that a system 
would have an equal number of binary symbols of each type, the average (over the long-term) 
data rate is nominally 13 kbps. If M-ary signaling is implemented with the same maximum and 
minimum symbol durations, the data rate can be improved by factor of log2(M), but at potentially 
increased channel BER. Achievable data rates are in principle limited by operating conditions 
(received SNR, tolerable BER, cumulant estimation errors, etc.).In addition to reliably detecting 
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a change in the SFOCMP eigenvalues using a basic correlation technique, a minimum b 
consecutive blocks are required (currently b=4) for each of S vector samples from the receive 
array, thereby making the theoretical minimum symbol duration equal to bS(l/R). Similarly, the 
incremental time duration for the alphabet should be at least ES(1/R), where E is the number of 
blocks desired by the designer to provide a balance between adequate safety margin in the time 
duration decision process and required throughput rate. In theory, E can be as low as unity 
which would enhance the achievable signaling rate for a fixed alphabet size. However, this is 
likely not a practical choice since numerous errors can occur due to the receive block processing 
not being time aligned to transmitter symbol boundaries. 

[0049] The receiver 303 uses an A/-element (or port) receive array 327 and an RF 
processor 305 to obtain the transmitted signal. In order to capture the temporal character (i.e. the 
time duration modulation of the SFOCMP eigenvalues) of the covert signal, the array data is first 
sampled and digitized at some rate suitable for the application. Each array output is digitized 
simultaneously producing a vector observation in the vector digitizer and buffer 307. The array 
output data is buffered and subdivided into non-overlapping blocks in 307. Block-wise across 
signal samples (i.e. the vector observations) are then collected from an array at the intended 
receiver aperture and the cumulants are block estimated, the matrix pencil is formed, and the 
generalized eigenvalue decomposition (GEVD) is performed by the Blind Source Separation 
processor 309. The operation of the BSS requires the selection of a triplicate of time lags 
provided by the time lags selection device 311. The GEVD provides a set of eigenvalues 

and N eigenvectors vjf^^ , where k=l, 2, 3, . , N (i.e. assuming an A^-port array is used) for 

each block of data. The superscript b is used as a block coimter in the receiver. We assume 
there are Ms generalized eigenvalues representing the SFOCMP properties for each of the Ms 
signals in the field of view (FOV) of the receive array 327, M <N , The remaining N-Ms 
eigenvalues are of the indeterminate (i.e. 0/0 type). Thus when using a sequence of block 
estimates for the SFOCMP eigenvalues of the Ms consistent signals will be apparent. Further 
processing is required and performed in the commxmication receive processor 319 to assemble 
valid messages. A major part of this processing is to use spatial information available from the 
GEVD processing. The spatial information regarding the transmitter location and the message 
content are linked in a 1 : 1 fashion by the generalized eigenvalues produced by the processing in 
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309. Using the "side information" of the available spatial variable greatly eases message 
recovery since we assume that a transmitter spatial location will be "slowly" varying (i.e., 
changing at a rate much less than the symbol rate of the message) , hence a message can be 
reconstructed in part by looking for message symbols represented by eigenvalues and their 
durations, associated with a "consistent" location. The designer must ensure that the symbol 
duration alphabet has sufficient minimimi support and increments such that the practical time 
duration recovery issues where ambiguous results can be obtained do not adversely affect the 
system performance. 

[0050] As may be apparent to those of skill in the art, there may be some advantage to 
overlapping blocks of the data. However, the following discussion deals with non-overlapping 
blocks. On each block, the two fourth-order spatial cumulant matrices required to form the 
SFOCMP are formed using pre-selected delay triplets. The delays can be either pre-selected, or 
subjected to online modification using a programmed search routine (if necessary). This search 
routine might be necessary when certain conditions, such as repeated eigenvalues for different 
signals are encoxmtered However, provisions are made for signals whose eigenstructure match at 
the delays selected to be repressed at different delays to provide improved discrimination if 
desirable. After the matrix pencil is formed, the GEVD is computed. From the GEVD, the 
eigenvalues and eigenvectors are used to determine the signal environment over time block b. 
Subsequently, the eigenvectors are used to determine the signal steering vectors and then the 
eigenstructure is correlated block wise in the Blockwise Eigenvalue Correlator 313 to determine 
any changes in the signal environment. A change, such as symbol boundary, in the 
communication signal will alter its contribution to the signal environment eigenstructure, 
measured by the SFOCMP, in a detectable maimer. This means a "significant" movement in the 
complex plane of eigenvalues. As signal changes are detected, those signals are cued for storage 
in the signal history database 317. As part of the storage procedure, the spatial location of the 
signal is determined (i.e., either angle-of-arrival (AOA) or geolocation, whichever appHes given 
the specific application) by the AOA/Geolocation processor 315. Additionally, the steering 
vector can be recorded, which is useful when refined spatial information is imavailable and the 
relative motion of the transmitter and receiver is negligible. The eigenvalues no longer 
correlating with the present signal structure are also written to the database. The temporal 
support (i.e. durations) of the eigenvalues no longer correlating with the current signal structure 
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is measured and stored. All this data is formed and recorded in the signal history database 317 
along with other ancillary data that may be useful for signal post-processing applications such as 
data mining or covert message recovery. 

[0051] An important function of a tracker is the track initiation and deletion logic. An 
embodiment of the tracks uses a fixed distance and a fixed number of consecutive "good 
associations" for initiation and a single "no association" for a track deletion. A "good 
association" is any measurement that is "close enough" to track. A "no association" condition 
occurs when all the measurements are "Too far" from a particular track, the distance indicative 
of an good association can be set empirically or experimentally. Track initiation and track 
deletion strategies can also be used to adapt to various situation. A Kalman-like approach to 
association gates can be adapted as the nvmiber of observation for a track are accumulated. Such 
an approach also has the advantage of replacing fixed averaging of the measurements. 

[0052] The design also allows for multiple access for communications. Consider the 
case where multiple remote covert emitters are sending data. It is unlikely that they would have 
exactly the same fourth-order cumulant representation, even if they are using the same base 
waveform. This is because any deviation from nominal waveform implementation (e.g., 
firequency change, waveform change, matrix pencil eigenvalue change , phase noise, I/Q 
imbalance, timing jitter, phase jitter, symbol rate change, pulse shape change, a fourth-order 
statistic change, relative rotational alignment of a signal constellation change, power amplifier 
rise/fall time change, and Doppler shift change) causes the fourth-order statistics of these signals 
to differ. Further, the multiple access signals are assumed distinguishable by spatial location. Of 
course this requires enough data to be collected to resolve the location, and the array must also 
provide such resolving power. But, if automated location is not possible at the receiver, due to, 
for instance, no calibration, the covert transmitters may still have multiple access if the multiple 
access signals can be assured uniqueness amongst themselves and the environment of sufficient 
degree in the SFOCMP eigenvalues. The receiver need not know the exact eigenvalues that will 
be used, but in this mode it is incumbent on the individual transmitters to use one and only one 
eigenvalue and not switch waveforms. In principle, correlation algorithms to properly sort this 
data are readily imaginable, though the details depend specifically on the signal designs. 
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[0053] Figure Sa is a spectrogram of the emitter message shown in Figure 4 when using 
the frequency shifts of the preferred embodiment discussed herein. The symbol boundaries 
defined by frequency offsets at the transmitter are clearly evident in the spectrogram. Figure 5b 
is a representation of a typical Electronic Support Measures (ESM) receiver/Signal Intelligence 
Spectral Monitor output viewing the received signal of Figure 5a. With reference to the 
spectrogram shown in Figure 5b, the time varying spectral pattern of the message shown in 5a is 
buried in interference and noise (SNIR<-10db) making the covert message very difficult to 
detect, and thus even more difficult to intercept and exploit. As is know to those of skill in the 
art the acronym SNIR stands for Signal to Noise and Interference Ratio. 

[0054] Given an environment with several interferers and the afready negative received 
SNR an unintended receiver (even using a front-end filter) will likely not reliably detect the 
presence of the covert signal. But even if a machine detects the presence of the signal energy, it 
would likely not be acted upon since it would fail all modulation recognition tests and show no 
exploitable temporal structure. The signal represented in Figure 4b is frequency shifted 
Laplacian (double-exponential) Noise. Viewed by a casual observer the signal would mimic 
additional thermal noise, hence even if detected by an energy detector of suitable design, the 
detection would likely be discarded. 

[0055] A mathematical element of the invention is the use of spatial high order statistics 
to separate signal sources, such as a blind source separation algorithm that utilizes a normalized 
spatial fourth-order cumulant matrix pencil and its generalized eigenvalue decomposition 
(GEVD). The equations presented herein use the following subscripting convention. Quantities 
relating to the array observations available to the system are denoted with a boldface subscript x. 
However, the subscript should not be confused with the representation of the vector observation 
from the array output, also denoted as a boldface x. From the context the meanings shall be clear 
to those of skill in the art. Further, quantities relating to the propagating signals impinging on a 
receive array are denoted with a boldface subscript r. Following this convention, the matrix 
pencil of the array output data is given as is given as equation 1. An assiunption is made that the 
received signals r comprising the vector observation of the array output x are independent. 
Therefore the spatial fourth-order cumulant matrix pencil (SFOCMP) of the array output can 
be written as: 
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10056] P, (A,t) = (0,0,0) - /ICj (ri,z-2,T3) (1) 

[0057] where the arguments of the pencil represent a generalized eigenvalue, X, and a 
triplet of time delays, x. The theoretical set of finite generalized eigenvalues turns out to be the 
inverse of the normalized fourth-order autocumulants of the M signals, [r^ (0}/=i the field of 

view (FOV) during the observation interval. The terms represent the spatial fourth-order 

autocumulant matrices. The arguments of the terms indicate the triplet of time delays used to 
form the matrices. The explicit computation is given as 




[0058] where the matrix isNxN, and the subscript rc indicates the element in the r 
row and the c^^ column. The subscript on the function x in the argument on the right-hand side, 
indicates which array port, i,r,c'=l,2.,.,Ny is being used. 

Because of the unique definition of the pencil of the array output data is related to the pencil 
of the impinging signals Pr as given in equation 2: 



P, (A,T) = (0,0,0) - ZCl{Ti,T2,T,) 

[00591 = v[c; (0,0,0) - ;iC;(ri,r2,r3)]v^ (2) 

= VPr(/l,T)V'' 



[0060] The quantity V shown in equation 2 is a ATx Af, matrix composed of the steering 
vectors for each signal impinging on the array, where N is the number of array ports available to 



WSH\96179.1 



13 



GCSD-1448 



HAR62.013 



the user and Af,, Af , < , is the number of signals. In a very simplistic and idealized case the 
well-known array propagation vector is a steering vector (i.e., the time delay is represented as 
phase). However, in general if the array is well-designed (i.e., no grating lobes) and the signals 
are emitted from non-identical locations, then the matrix V is of full rank. This guarantees an 
equivalence between the eigen structure of the pencils Pr and P,. 

[0061] Since Pr is a pencil solely of the received signals, and the signals are assxuned 
independent, then by virtue of the properties of cumulants, the pencil Pr is diagonal. This 
property does not hold true for the pencil formed with the array output data x. However, because 
of "equivalence'' finite eigenvalues of Px are the finite eigenvalues of Pr, access to an exploitable 
high-order statistical property the eigenstructwe of the SFOCMP is available. As introduced 
here these eigenvalues represent the fourth-order characteristics of each received signal. . 

Specifically, each signal in {r,- (0}/^i contributes one finite eigenvalue, and it is expressed as the 
inverse normalized fourth-order autocumulant for that signal as expressed by equation 3. 



ct (0,0,0) 

[0062] ^=-7^ for/w = l,2,...,M (3) 



4 

[0063] where the terms represent the individual fourth-order cumulant terms for each 

m 

signal. These terms are actually the diagonal terms of the pencil Pr as shown in equation (4). 



WSH\96179.1 



14 



GCSD-1448 



HAR62-013 







« c • • • • 




0 


« 




• 
• 


c;(oiQO)-;^(Ti.^,t^) 




0 





(4) 

[0064] Thus the GEVD of the two pencils and Pr have the same set of finite solutions 
for the eigenvalues. The eigenvalues are the terms where the rank of the pencil is reduced. It 
should be readily apparent that values given by equation (3) are the eigenvalues of the pencil 
equation (1). 

[0065] These eigenvalues are available to an analysis system, and in theory are 
independent of system Gaussian noise level given sufficient length data records. The 
eigenvalues are implicit characteristics of the signals carrying the emitter's covert message in 
each symbol duration. To exploit this property, as mentioned before, the receiver will typically 
fomfi blocks or batches of received data for the purpose of correlating the eigenstructure over 
time to determine patterns of persistent values (Figure 3) augmented by the availability of spatial 
data. It is important to recall that only the time duration of the emitter's statistical characteristic 
as measured by the SFOCMP is relevant, and not the exact values. Hence, the emitter is 
completely free to choose the "carrier waveforms" at will (Figures 6a-f). 

[0066] The steering vectors can be estimated from the cumulant data for each signal in 
the FOV of the receiver. A cumulant matrix formed by the receive data, say (0,0,0) and for 
each eigenvector available from the pencil P, forms. 



[00671 C^(0,0,0)e^/^ = 



i;c,^(0,0,0)v,v/ 
.1=1 
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[0068] The last equality follows directly from the fact that each eigenvector of the 
SFOCMP Px is orthogonal to each signals steering vector, v ."e^/^ = 0 when j , This fact is 
generated by the unique construction of the SFOCMP and the definitions of the cumulants. 

[0069] In Figm-e 3, the Blind Source Separation processor 309 forms and applies a 
separation Matrix and enumerates the number of sources. As described above, from an array 
output the spatial 4^^ order cumulant matrices are estimated and the estimates are used to 
determine the eigen analysis for the first-order matrix pencil. Signal detection and enumeration 
providing the number of source is performed and the separation matrix from the pencil 
eigenvectors is accomplish. The separation matrix is also used to determine spatial variables, 
such as steering vectors, AoA or geolocation in block 315 of Figure 3. 

[0070] There are numerous ways the covert transmitter can control the desired 
characteristic of its emitted waveform, some of which leave lower-order statistics unchanged. 
For example altering the channel filter (i.e. Nyquist pulse shaping) between maximum phase and 
minimum phase realizations is undetectable in the second-order domain (i.e. power spectra), but 
evident in the fourth-order domain as measured by our SFOCMP. Also one could conceive on 
signaling with kurtosis, a fourth-order statistic, applied in the transmitter. Or, one could simply 
shift between classic waveforms, for example, BPSK, OMSK, QPSK, QAM, or potentially even 
just variants (i.e. constellation rotations, different pulse shape filters) of a fixed modulation type. 
There also exists the possibility that the "carrier waveforms" might be chosen as chaotic to 
appear more noise-like or designed using numerical techniques and generated using direct 
synthesis in a transmitter. 

[0071] In any scheme adopted, the information transfer from the transmitter to the 
receiver is contained in the duration of the change in the eigenvalues of the SFOCMP, not the 
particular eigenvalues. Since our technique is independent of the particular eigenvalue, it is 
independent of the waveforms used by the emitter. Which allows, in principle, any transmitter to 
make use of the receiver having the capability to exploit fourth-order cumulants. The degree to 
which a specific emitter wishes to "hide" from say conventional ESM receivers holds the 
implications on the implementation details of the "carrier waveforms". 
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[0072] Figure 6a is a generalized schematic diagram of an embodiment of a noise signal 
generator for a waveform independent covert conmiunication system. A noise generator 602 
generates temporal dependent non-Gaussian noise. The output of the noise generator 602 is 
combined by combiner 606 with a carrier signal or waveform source 604. The combiner can be 
an adder, a mixer, a multiplier, a non-linear device or other type of combiner that facilitate a 
change in a higher order signal statistic. This Ust of combiners in not exhaustive and should not 
be construed to limit the scope of the invention. The combined signal forms a type of baseband 
signal that is processed (eg. amplified, upconverted, etc.) and transmitted by transmitter 608. 
The generalized elements in Figure 6a form a basic framework and facilitates several different 
embodiment for noise signal generation and transmission of a covert message. The generalized 
diagram in 6a as well as embodiment shown in Figures 6a-6e presented to illustrate variations of 
noise signal generators and not intended to limit the scope of the disclosure. 

[0073] Figure 6b is a specific implementation of Figure 6a with dual signal sources. The 
noise generator 602 is implemented with a Non-Gaussian Noise Generator (NGNG) 610 and a 
temporal filter 612 that creates an output with temporal dependence. The signal source 604 is 
implemented with two unique signal sources 614a and 614b. Unique signal sources 614a and 
614b are selected to create unique spatial high order statistics when combined by the combiner 
606 with the output of the noise generator 602. The unique signal sources 614a and 614b are 
connected to the combiner by a switch 618. The switch is driven by the symbol bit stream of the 
covert message. The switch 618 connects altemate unique signal sources at the conclusion of 
each successive symbol duration. It is important to recall that only the time duration of the 
emitter's statistical characteristic as measured by the SFOCMP is relevant , and not the exact 
values. Therefore the emitter is completely free to choose the carrier waveforms. While two 
unique signal sources are shown in Figure 6b, multiple unique signal source can likewise be 
applied with the switch altemating between the signal sources. 

[0074] Figure 6c is another specific implementation of Figure 6a having dual temporal 
filters. The noise generator 602 is implement with a non-Gaussian noise generator 610 
connected to two unique temporal filters 612a and 612b. The unique temporal filters 612a and 
612b are selected to create unique spatial high order statistics of the transmitted combined signal 
(e.g. maximum and minimum phase filters). The unique temporal filters 612a and 612b are 
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connected by switch 618 to the combiner 606 and combined with the output of the signal source 
604. Again the switch is driven by the symbol bit stream of the covert message in the same 
manner described above for the dual signal source implementation. The signal or waveform 
transmitted from transmitter 608 of Figure 6c provides altemating unique spatial high order 
statistics for each successive message symbol independent of the symbol transmitter. 

[0075] Figure 6d is a specific implementation of the noise generator 602 including two 
unique non Gaussian noise generators 610a and 610b connected by switch 618 to the temporal 
filter 612. Again the output of the noise generator combined with the signal soiu^ce 604 in 
combiner 606 and transmitted by emitter or transmitter 608. The imique NGNGs 610a and 610b 
are likewise selected for their effect on spatial high order statistic of the transmitter signal. 

[0076] Figure 6e is an implementation similar to the noise generator 602 of Figure 6d, 
where a non-Gaussian noise generator 610a and temporal filter 612a form a imique branch and 
NGNG 610b and temporal filter 612b form another unique branch. The branches are unique in 
the sense of their respective spatial high order statistics. The branches are connected to the 
combiner 606 by switch 618. 

[0077] Combinations of the specific implementation described and others that should be 
readily apparent from an understanding of this disclosure are likewise envisioned. Figure 6F is 
but one of the many possible combinations. The implementation of the noise generator 602 of 
Figure 6f includes the noise generator implementation of Figure 6d containing two unique 
NGNGs 610a and 610b coupled with the signal source 604 of Figure 6b containing two unique 
signal sources 614a and 614b. While the switches 618a and 618b are constrained to switching at 
only at the duration boundaries, their operation advantageously would operate independently. 
The result of this specific implementation allows for switching between four unique signals. 

[0078] Figure 7 is an embodiment using a Laplacian generator 710 and an infinite 
impulse response low pass filter 712 (e.g. an autoregressive single-pole filter) as the noise 
generator 602. The Laplacian generator 710 is chosen for leptokurtosis property. We have 
chosen to form the necessary sequence of random Laplacian variates by selecting a pair of 
random variants, say x and y, uniformly distributed between 0 and 1 and forming log(x/y) for 
each sample required. For complex baseband symbols a total of four uniformly distributed 
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variates must be sampled providing a stream of complex samples where the real and imaginary 
parts of each are independently Laplacian distributed. The output of the noise generator 602 is 
combined by combiner 606 with the signal source 604. The signal source 604 is implemented 
with multiple unique signal sources 714i, 7142,. • .714i that are used to switch spatial high order 
statistics (i.e. matrix pencil eignevalues for the present embodiment) for each symbol of the 
covert messages. The unique signal sources 714 shown are offset in frequency to provide unique 
matrix pencil eigenvalues. In this particular embodiment, L must be greater than or equal to the 
alphabet size M. 

[0079] As shall be imderstood by those of skill in the art, the specific example discussed 
above may be extended to use random mappings of frequency offsets over time. Also, we could 
alter the channel filter. There is no requirement that the filter be IIR as shown in the Figure. A 
number of altemative implementations could be chosen depending on the application. The key 
feature of the filter is to introduce a temporal dependence of the input noise waveform. Further 
one could also consider altering the input noise generator. However, a consideration is to select 
a source with suitable fourth-order properties. Any or all of these parameters can be modified to 
control the fourth-order properties for the transmitted waveform so long as the "codebook" 
constraints (time duration and alphabet size) are maintained. A natural altemative to frequency 
shifting, would be to pulse the carrier on/off. However, this approach reduces the nximber of 
signal samples available for geolocation given a fixed observation interval as discussed hereafter. 

[0080] If one wished to use standard waveforms as the "carrier waveforms" this mode of 
operation is also possible with this invention. The transmitter shown in Figure 7 could be 
modified to look like alternatives shown in Figures 6b-6f It is also important to notice that any 
deviations from a nominal waveform type, such as a QPSK waveform without phase noise or 
timing or I/Q imbalance, will cause a detectable shift in a signal's fourth-order statistic. In 
principle, the basic waveform need not be altered from QPSK to say BPSK, but uniqueness can 
be achieve by controlling the parameters of, for example, timing jitter, "carrier waveform" 
symbol rate, phase noise, pulse shaping and the like to provide the necessary separation in the 
fourth-order domain or other high order domains. The "carrier waveforms" can contain 
completely worthless symbols, hence anyone trying to conventionally exploit the "carrier 
waveform" information would be wasting their time, alternatively the "carrier waveforms" may 
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carry information with a secondary message. However, a secondary message contained within 
the carrier waveform would not necessarily be LPE, LPI or LPD. 

[0081] An example of a potential message recovery embodiment is shown in Figure 8. 
To recover a message the receiver requires three parameters relative to the time durations of the 
message symbols. The receiver must know the minimum symbol duration 801, maximum 
symbol duration 803, and some other information regarding the other M-2 symbols in the 
alphabet. In the embodiment shown in Figure 8 only a binary symbol alphabet is used. This 
information can be provided, say by policy where the symbol duration increments are integer- 
related to the minimimi symbol duration 801. In this case the receiver need only know the 
maximum duration, the minimum duration and alphabet sizes. Otherwise the receiver can just 
know a-priori the pre-selected durations which can then be arbitrarily selected. The symbol 
durations should be selected in an application for easing the decision process in the receiver to 
map the measured durations to one of the alphabet elements. In this case we used a symbol 
twice as long as another, but such widely disparate time durations may impact throughput. A 
minimum symbol duration threshold 807 and a maximum symbol duration threshold 811 can be 
used as a time gate to remove spurious or interfering tracks. Spatially correlated signal tracks 
823, 825, 827 and 829 are examples of eigenvalue durations which exceed the threshold in the 
time gate and thus are ignored or discarded. 

[0082] Spatially correlated eigenvalue time durations 8059, 805 16, 80536, 8O553, 80568 and 
80580 are sequenced in signal track 821. The eigenvalue time duration of signal track 821 are 
compared to a decision threshold 809 to map and recover the encoded message "101011" of 
signal track 821. Where a M-ary alphabet is used M-l decision thresholds are required. 

[0083] An illustration of a portion of the block-to-block eigenvalue correlation result is 
shown in Figure 9. For block 30 of the covert emitter example in Figure 4, the large complex 
plane diagram in the top portion of the Figure 9 shows the complex eigenvalue locations of the 
SFOCMP (GEVD) results and the predicted locations of the blockwise eigenvalue correlator. 
The legend identifies the four levels of eigenvalue correlation confidence, ('"New," "Tentative," 
"Candidate," and "Confirmed"). The five consistent signal eigenvalues of the four co chaimel 
interferers and of the covert emitter are indicated by the smaller rectangle. The inconsistent non- 
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signal eigenvalue outside this box tends to move about the complex plane in an erratic/unstable 
fashion from one block to the next. In this block 30 of this example, the covert emitter 
eigenvalue lies below the imaginary axis. Blocks 30 and 31 correspond to two different 
symbols. The lower portion of Figure 9 illustrates the blockwise changes in eigenvalue locations 
over blocks 30 to 34 where the message symbol "boundary" occurs between blocks 30 and 31. 
For blocks 31 to 34, the covert emitter eigenvalue Ues in a consistent location above the 
imaginary axis of the complex plane. The stronger co channel interference source eigenvalues, 
within the oval for block 34, Ue in consistent locations over all five blocks. In a case involving 
pulsed interference emitters, the consistency of the spatial locations of the interference emitters 
and covert emitter allow message recovery embodiment to "stitch" together eigenvalues that fall 
within a "time-gate" and come from the same spatial location. 

[0084] As mentioned above, using a simple time-gating operation in the receiver, it is 
possible to determine which eigenvalues are potentially information carrying. By correlating the 
GEVD over many blocks of data the persistence of the eigenvalues can be measured. The 
persistence of eigenvalues of the SFOCMP over time from the covert transmitter provides the 
signaling mechanism. However, there may be a number of extraneous pulsed signals in the FOV 
time coincident with the desired communication signal This makes message recovery 
complicated, though with proper message construction and error recovery/correction (i.e. FEC), 
the system is robust to several types of errors such as "erasures" when ambiguous results may be 
obtained in the decoding and symbol recovery errors. These results can be encoxmtered due to 
signal fades (i.e. erasures) or symbol recovery errors in the receiver due to statistical fluctuations 
in time duration measurements exceeding a tolerable threshold. We can correct improper 
decisions regarding the detection of a symbol in the message in the receiver using typical error 
control coding. 

[0085] An embodiment of the receiver 303 is shown in Figure 10. The receiver 303 
augments the message recovery process by using spatial information regarding the locations of 
the transmitters in the field of view. In Figure 10 signal detections are performed in processor 
1001 which is composed of a BSS processor 309 and (eigenvalue) tracking processor 313 
previously described. The results of the blockwise processing and eigenvalue correlation 
previously described are stored in the Signal History Database 317. The data envisioned for the 
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database are signal characteristics measured which include but are not limited to eigenvalue 
identification number, time initialized, time deleted, number of observations, an estimate of the 
eigenvalue, and importantly an estimate of the spatial location associated with the eigenvalue. 
To recover the message in the communication receive processor 319, it is advantageous to 
initially time gate (filter) all the detections (i.e. durations of the eigenvalues) to limit the scope of 
the message recovery process. The detections surviving the time gate test are stored in the signal 
history database 317 along with other associated eigenvalue characteristics, including but not 
limited to Signal track identification. Symbol estimate. Time initialized, time deleted, number of 
observations and importantly spatial location. The detections are assembled into strings of 
symbols (tracks) using the criteria that a message must associate with a "consistent" spatial 
location. The preferred spatial information is angle-of-arrival (AoA) or geolocation for each 
signal detected in each block which are determined in the AoA/Geo computation processor 315. 
This requires that the receiver platform have knowledge of the iV-port sensor array calibration 
and the platform position during receipt of the covert message, such knowledge would normally 
be available. For simplicity of presentation, firom hereon the term '^spatial location" or spatial 
variable refers to AoA if only calibration is known, or geolocation if both receiver position and 
array calibration are known. It is preferable that the spatial relationship between the transmitter 
302 and receive array 327 be fixed over the message length, however, a slowly varying spatial 
property can be accommodated by introducing a method to correlate a sequence of spatial 
variables across contiguous blocks of data in the receiver 303 specifically in the Message 
Sequence Correlator 323 of the Covert Conmiunication Processor 319. 

[0086] Spatial correlation can be broaden to include simply steering vectors. This is 
usefiil when the array and transmitter have a stable geometry. Relative motion between the 
transmitter and sensing array causes the steering vectors to have a detrimental time dependency. 
Again, if the spatial variable for correlating the message data is "slowly" varying then small 
incremental changes can be tolerated. The covert messages is indented to be recovered using a 
"consistency" of the spatial domain information of the computed eigenstructure fi'om the 
SFOCMP for the signals of interest. But to account for the possibility of "fixed" location 
emitters and other emitter who are not of interest to the conmiunication process, a time-gate 
decision process as noted earlier is advantageously applied. This way the receiver need only 
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attempt to decode "message strings" that emanate from "consistent" spatial locations with the 
appropriate time character. 

[0087] Although access to the spatial variables using only the receive array output data 
has been previously described. It is useful to note a blind source separation algorithm based on a 
fourth-order cxmiulant matrix pencil produces eigenvectors that are orthogonal all but one 
signal's steering vector. Thus using the eigenvectors it is possible to estimate each 
corresponding signal's steering vector. Two methods are possible. The first method is to use the 
blockwise estimates directly available from the BSS process as described in relation to Figure 3 
such as averaging across the blocks. A second method is to use the time stamps available from 
the blockwise correlation process 313 of particular eigenvalues as indexes into the digitizer 
buffer memory. Subsequently, the raw data so indexed is used to directly make an estimate of 
the steering vector. In both cases the steering vector will be estimated using the data available 
from the GEVD of the SFOCMP. Once steering vector estimates are available, the estimation of 
the other spatial variables, AoA and/or geolocation, can be determine by methods well-known in 
the art. The characteristic that allows this computation is that the eigenvalues and eigenvectors 
of the GEVD of the SFOCMP have a 1 : 1 correspondence as in standard eigenanalysis. So, when 
the eigenvalues are used to measure the high-order statistical properties of the received signals, 
an index relating directly to where that signal with that characteristic emanated from is 
available. Again, the spatial dimension for signal association can be exploited to remove any 
ambiguity of the temporal decomposition, since we assume that no two emitters are identically 
located. 

[0088] Signal tracks 1021, 1023 and 1025 containing associating matrix pencil 
eigenvalues from consistent sources as evident from common spatial variables are used to 
recover the message by correlating track information against spatial variables and time gates. 
This recovery process in performed by the covert conraiunication receive processor 319. 

[0089] The recovered messages are composed of sequentially ordered matrix pencil 
eigenvalues with a duration within the time gate originating from the same location as 
determined by the steering vector estimate, AoA or geolocation. The "signal sequence 1" 1021 
includes eigenvalues 1, 4, 7, 10, and 12 with durations mapping message "0101 1", these 
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eigenvalues and their associated signals all originated from a "consistent" spatial location, 
namely, Geoi, AoAi or had the same steering vector SVi. Similarly, "signal sequence 2" 1023 
includes eigenvalues and durations with common steering vectors, AoA, or geolocation and the 
durations pass the time gate filter thus recovering a "000100" message. "Signal sequence 3" 
1025 however had two eigenvalues 3 and 9 with high spatial correlation, however the signal is 
not decoded since the durations exceed the time gate allowable. In short, the receive processor 
recognized that the time durations in sequence 1025 do not correspond to any symbols of the 
expected alphabet. 

[0090] The advantages of incorporating spatial variables into the message recovery 
process warrant explanation. First, the spatial variables aid in rejecting extraneous pulsed 
emitters based on their spatial locations being anti-correlated over time to the persistent spatial 
locations of the covert emitter(s). By the same token, spatial variables allow the basic signaling 
approach to support multiple access of covert emitters without undue burden in the receiver for 
properly assembling the pulsed message sequences. This is because the additional covert 
emitters will very likely emanate from resolvable spatial locations, and the receiver can use the 
consistency of the spatial locations over time to associate the proper message sequence. For each 
transmitted signal, the message sequence is represented by the time durations of the eigenvalues 
of the appropriately selected matrix pencil, where we have preferred the SFOCMP approach. 

[0091] The spatial location of any emitter is independent of the exact value of its 
corresponding eigenvalues available from the GEVD of the SFOCMP. Lastly, the spatial 
variables provide additional "distance" in the recovery process, since it is now multi- 
dimensional. For example, two signals may have very similar eigenvalues. But, if their spatial 
locations are resolvable by the receiver, and fairly constant, then the eigenvalues corresponding 
to those spatial locations can be easily assigned. Then the message can be recovered using the 
time duration of each eigenvalue in the sequence assigned a given spatial location using the same 
technique as previously described when only a single covert signal was in view. 

[0092] The ability to resolve spatial location has system implications that are interrelated. 
Some top-level practical design issues that must be reconciled are desired proximity of 
transmitters, expected noise environment, block processing sample support for estimating spatial 
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locations and the eigenvalues available from the SFOCMP, digitizer sample rates, signal 
bandwidths and center frequencies, aperture design (i.e. element type, size, number of ports, 
operating frequency), and the like. This is of course in addition to having an appropriate level of 
calibration and positional knowledge of the receive platform. Many of these considerations are 
direct carry-overs from standard array-based signal processing systems. 

[0093] Figure 1 1 is a flow chart for covert communication with a transmitter 1 102 and a 
receiver 1 103 using message recovery with spatial information as shown in Figure 10. The 
source 1 100 provides a data stream to the covert transmitter 1 102 as a stream of M-ary alphabet 
symbols. In block 1 104 each M-ary alphabet symbol is assigned a unique duration. Plural 
waveforms with determinable high order statistic that are constant and non zero are generated in 
block 1 106. One of the plural waveforms is selected in block 1 108 and is transmitted for the 
symbol assigned duration in block 1110. In block 1 1 12 it is determined if the transmitted 
symbol is the last in the message, if the last symbol has be transmitted the process is returned in 
block 1 1 16 to standby or other quasi active state pending initiation of a new message. If other 
symbols remain to be transmitted, a different waveform from that previously transmitted is 
selected in block 1114 and transmitted for the respective symbols assigned duration in block 
1110. 

[0094] The receiver 1 103 receives the transmitted waveform along with environmental 
and random noise in a multi element array as shown in block 1 105. Using BSS and GEVD, 
matrix pencil eigenvalues (MPE) are obtained in block 1 107 and the MPE are blockwise 
correlated to determine their respective durations in block 11 09. A spatial variable or location is 
determined from the MPE in block 1111, again the spatial variable can be selected from a 
Steering vector, AoA or geolocation. The MPE detection durations are time gate filtered in 
block 1 1 13 to reject detections outside of the duration minimum and maximum thresholds. The 
passed MPE and durations are sequentially sorted by common or consistent spatial variable in to 
signal tracks in block 1115. The signal tracks are maps are mapped ( i.e. their durations are 
correlated to symbols and their respective assigned durations)to recover the message in block 
1117. The message may be processed otherwise transformed to provide the message at the sink 
1101. 
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[0095] If spatial data is unavailable, say because calibration of the sensing array has been 
degraded, the communication process can still operate. However, the freedom of waveform 
selection by the transmitter is reduced. In this case the transmitter must select a specific 
waveform type and use it exclusively (in a pulsed fashion) over the entire message (Figure 4). 
Unfortunately with this implementation option the achievable data rate is reduced because of the 
need to introduce "deadtime" to define symbol boimdaries. Thus it is possible to use a single 
"carrier waveform" that can be pulsed "ON" for each symbol followed by a period of "OFF" 
time. In this way the covert transmitter need only use a single waveform and need not modify 
it's fourth-order cumulant signature. This could be an advantage in systems where additional 
spatial correlation variables preferred to aid unambiguous assignment of the received 
eigenvalues are unavailable, since we identify one and only one eigenvalue. But as mentioned 
this is disadvantageous in a multi-emitter environment. For such an environment we prefer a 
waveform agile emitter where the "carrier waveform" sequencer logic can be designed to select 
the "carrier waveform" for a specific duration, controlled by the message symbol, in either a 
fixed map or in some other manner. The mapping choice would be up to the transmitter designer 
and need not be known to the receiver. The receiver processing is shown in Figure 12. 

[0096] Figure 12 presents a receiver 303 for processing covert messages without spatial 
location or spatial variable information. Processor 1201 contains BSS processing 309 and 
eigenvalue correlation 313 operating in the same manner as previously described, however 
without transferring or developing spatial variables. Therefore the eigenvalue characteristic and 
time stamps for determining eigenvalue duration are stored in the signal history database 317, 
without a spatial variable. For the embodiment of Figure 12, the tracks are not sorted by spatial 
location by rather eigenvalues. Signal tracks contain only detections of the same eigenvalue with 
a detected duration that passes the time gate. Signal 1221 includes only one eigenvalue with 4 
detection durations that pass the time gate. Signal 1221 thus is mapped to recover message 
"0010". The durations of the "OFF" period of pulse is irrelevant as long as its duration is greater 
that 1 block to enable the previous detection to be terminated the tracking correlator 313. Signal 
tracks 1223 and 122S are both ignored or discarded for failing the time gate. 

[0097] Figure 13 is a flow chart for covert communication with a transmitter 1302 and a 
receiver 1303 using message recovery without spatial information as shown in Figure 12. The 
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process is similar to that shown in Figure 1 1 . The source 1300 provides a data stream to the 
covert transmitter 1302 as a M-ary alphabet symbols. In block 1304 each M-ary alphabet 
symbol is assigned a unique duration. However, now a signal waveform with a determinable 
high order statistic that is constant and nonzero is generated in block 1306. The waveform is 
transmitted for the symbol assigned duration in block 1310. In block 1312 it is determined if the 
transmitted symbol is the last in the message, if the last symbol has been transmitted the process 
is returned in block 1316 to standby or other quasi active state pending initiation of a new 
message. If other symbols remain to be transmitted, the transmitter is paused, or another imique 
waveform is transmitter for an small duration t=k in block 1320. The waveform generated in 
block 1306 is again transmitted for the respective symbol assigned duration in block 1310. 

[0098] The receiver 1303 receives the transmitted waveform along with environmental 
and random noise in a multi element array as shown in block 1305. Using BSS and GEVD, 
matrix pencil eigenvalues are obtained in block 1307 and the MPE are tracked to determine their 
respective durations in block 1309. The MPE detection durations are time gate filtered in block 
1313 to reject detections outside of the duration minimum and maximxmi thresholds. The passed 
MPE and durations are sequentially sorted by specific eigenvalue into signal tracks in block 
1321. The eigenvalue durations forming the "signal sequences" are mapped (i.e., their durations 
are correlated to symbols and their respective assigned durations) to recover the message in 
block 1317 and delivered to the sink 1301. 

[0099] Figure 14 is a flow diagram for covert communication using encoded eigenvalues 
and spatial information. The source 1400 provides a data stream to the covert transmitter 1402 
as a M-ary alphabet symbols. In block 1404 the M-ary symbols of the messages are encoded as 
MPE with unique durations. The eigenvalue is transmitted for the symbol assigned duration as a 
respective source waveform in block 1410. The source waveform is the low order signal whose 
high order statistic creates the respective eigenvalue. In block 1412 it is determined if the 
transmitted symbol is the last in the message, if the last symbol has be transmitted the process is 
returned in block 1416 to standby or other quasi active state pending initiation of a new message. 
If other symbols remain to be transmitted, a different MPE is selected in block 1414 and 
transmitted with a respective symbol assigned duration as a respective source waveform in block 
1410. 
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[00100] The receiver 1403 receives the transmitted waveform along with environmental 
and random noise in a multi element array as shown in block 1405. Using BSS and GEVD, the 
matrix pencil eigenvalues encoding the symbol are recovered and tracked in block 1407 to 
determine their respective durations and spatial information in block 1409. The MPEs are 
decoded based on MPED and Spatial information as previously described in block 1417 to 
recover the message and provide the message symbols to the sink in block 1401. Generally the 
noise is white Gaussian noise, color noise or interferer signals. 

[00101] Figure 15 is a representation of a binary message sequence 101011 transmitted via 
data carrying waveforms according to an embodiment. In the embodiment shown the symbol are 
transmitted by alternating waveforms of BPSK, QPSK and GMSK. As the message is 
independent of the encoded data, for LPI it is preferred that the carrier waveforms be modulated 
with random data at a rate greatly exceeding the covert message symbol rate. The random data 
may be modulated for a BPSK waveform at 5 sample/per random symbol with a frequency offset 
of 0. 1 fs, for QPSK wave at 10 samples/random symbol with a frequency offset of -.05 fs and for 
a GMSK waveform (GSM:hmod=l/2, Btb=l/3) 2 samples/random symbol with a frequency offset 
of 0.0005 fs. Other waveforais including to DQPSK, DBPSK, FSK, QAM, DPCM are equally 
applicable and are also envisioned, however embodiments of the invention should not be 
construed to be limited to the particular waveforms listed. 

[00102] As discussed previously the waveform duration discriminates the message 
symbol. As shown in Figure 15 the symbol 1 has a duration of 20 blocks and the symbol "0" 
has a duration of 10 blocks. The message is transmitted as a GMSK waveform for a duration of 
20 blocks indicating a symbol "1". The next symbol is transmitted as a different waveform 
QPSK for a duration of 10 block indicating a symbol "0". The next symbol is transmitted by the 
GMSK waveform for a duration of 20 blocks again indicating a "1" symbol and a QPSK 
waveform is used to transmit the "0" symbol. As illustrative of the independent of the waveform 
and message content the symbol "1" is then transmitted as a BPSK waveform for a 20 block 
duration and the next symbol "1" is transmitted by the QPSK waveform for a duration of 20 
block. Evident in Figure 15 is that identical waveforms can communicate different symbols 
while maintaining waveform content and message independence. 
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[00103] The use of the higher-order statistics can be used to geolocate a transmitter. 
Figure 16 is a embodiment of a system for geolocation 1600, A receiver 1609 with a multi- 
element sensor array or antenna array receives the transmitted signal(s) 1621 from the target 
transmitter(s) 1624. The receiver 1609 contains processors, microprocessor and/or logic circuits 
implemented with hardware and/or software to estimate a spatial 4^^ order cummulant matrix 
1605, perform eigen analysis of the first order matrix pencil 1607, perform signal detection and 
enumeration 1608 and form the separation matrix 1611. The steering vectors are then estimated 
using the non-orthogonal eigenvectors with a 1:1 mapping of the eigenvalues obtained with the 
GEVD of the matrix pencil as shown in block 1620. Using calibration data of the multi-element 
sensor array the AoA of the detected signal(s) from the transmitter(s) in 1630. The AoA coupled 
with the position of the platform (multi-element sensor array) enables geolocation of the target 
transmitter(s) shown as block 1650. The determination of geolocation from AoA and positional 
data can be achieved from many known methods. 

[00104] The subject matter regarding geolocation an also be used to "track" a mobile 
convert emitter, or a mobile emitter with a distinct and fixed temporal structure or characteristic 
(i.e. a chain of eigenvalues) known to the receiver. The receiver uses GEVD indexed by the 
eigenvalues to compute and associate a sequence of AoAs and geolocations. A receiver with the 
a-priori knowledge of the temporal structure can associate a sequence eigenvectors with 
correlated temporal characteristic and applying kinematics constraints (i.e. maximum velocity 
and/or maximum acceleration for the emitter), construct a path history or prediction for the 
emitter. 

[00105] Figure 1 7 is a flow chart for tracking a cooperative or known emitter. The target 
emitter 1701 whether cooperative or non-cooperative transmits a signal with a known, fixed 
temporal characteristic in block 1709. The receiver 1702 receives the transmitted signal and 
other unknown signals collectively candidates signals with a multi-element sensor array, 
digitizes the output of the array for each of the candidate signals in block 1710. The blind source 
selection process discussed earlier is performed on the candidate signals, including determining a 
spatial 4^^ order cummulant 1711, determining a matrix pencil eigenvalue 1712 and performing 
generalized eigenvalue decomposition 1713. 
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[00106] From the generalized eigenvalue decomposition 1713 a spatial variable is 
determined in block 1715. A non-orthogonal eigenvector corresponding to the steering vector of 
the candidate signal is selected as the spatial variable. The receiver also determines a temporal 
characteristic of the candidate signals in block 1714. The spatial variable and the temporal 
characteristics of the candidate signal are associated in block 1716 and the temporal 
characteristics are correlated, or compared to the unique known temporal characteristic of the 
signal transmitted from the target transmitter in block 1717. The spatial variable associated with 
the temporal characteristics that are highly correlated with the temporal characteristics of the 
target transmitters signal are selected in block 1718. These selected spatial variables like AoA or 
geolocation can be used to track the target transmitter. 

[00107] While preferred embodiments of the present inventive system and method have 
been described, it is to be understood that the embodiments described are illustrative only and 
that the scope of the embodiments of the present inventive system and method is to be defined 
solely by the appended claims when accorded a full range of equivalence, many variations and 
modifications naturally occurring to those of skill in the art from a perusal hereof 
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